%P1误差
d=[0.5,0.25,0.125,0.0625,0.05,0.03125,0.02]
error1=[0.0241648,0.00671632,0.00189399,0.000473413,0.000304804,0.000118138,4.83182e-05]
%plot(log(d),log(error1))



%P2误差
d=[0.5,0.25,0.125,0.0625,0.05,0.03125,0.02]
error1=[0.000325077,4.49184e-05,5.83402e-06,7.16864e-07,3.57322e-07,8.59558e-08,2.20749e-08]
plot(log(d),log(error1))
